Waveform sample method of excitable sensory neuron 
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We present a new interpretation for encoding information of the period of input signals into spike- 
trains in individual sensory neuronal systems. The spike-train could be described as the waveform 
sample of the input signal which locks sample points to wave crests with randomness. Based on 
simulations of the Hodgkin-Huxley (HH) neuron responding to periodic inputs, we demonstrate that 
the random sampling is a proper encoding method in medium frequency region since power spectra 
of the reconstructed spike-trains are identical to that of neural signals. 
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In the last decade the phenomenon of the stochas- 
tic resonance of excitable neuron has been extensively 
studied both experimentally Q, 0, 0, H and theoretically 
H, IE I n these studies stochastic excitable sensory 
neurons or neuron models subject to periodic signals. 
Stochastic resonance means that the ability of neurons 
for responding to the weak signal is enhanced by noise. 
The power spectral density of action potential was used 
as one of common methods for estimating the sensitiv- 
ity of stochastic neurons to weak input signal 8]. The 
power spectrum revealed the frequency information of 
input signal from the time series action potential 0, |t| . 
However, the code for encodin g an d decoding in neuron 
is not clarified at the moment [lOj . A approximate the- 
ory for modelling neuron firing in the presence of noise 
and a periodic stimulus has been proposed and the 
process of neuron conveying information was described 
by inhomogeneous Poisson point process. In this paper, 
we didn't treat the spike-train as a poisson process. Our 
source of inspiration is based on the phase locking [Tlj . 
We investigated the correlation exhibited in time axis 
between the temporal sequences of firing events and the 
input signals, and proposed the waveform sample method 
used by neurons for processing periodic signals. 



three for the gating variables: m, n, and h; that is, 

V = (l lon (t) + Iext(t) + Inoise)/C m , (la) 

rh={m ao {V)-m)/T m (V), (lb) 
h={h oc (V)-h)/T h (V), (lc) 
h=(h oc (V)-n)/T n (V), (Id) 

where 

hon{t) = -9nam 3 h{V-V na )-g k n 4 (V-V k )-gi(V-Vi). 

(le) 

The ionic current Ii on (t) describes the ionic transport 
through the membrane, and includes the sodium (I na ), 
potassium (Ik), and leak (/;) currents. The constants 
9na, 9k, and gi are the maximal conductances for ion and 
leakage channels, and V na , V k , Vi are the corresponding 



reversal potentials; m c 



and T m , T h , r n rep- 



resent the saturation values and relaxation times of the 
gating variable s [Till . Detailed values of parameters can 
be found in [l2lll5|. L T> is the external stimulus received 
by neuron. We model the noisy current as an exponen- 
tially correlated colored noise 
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We made a study based on a popular Hodgkin-Huxley 
(HH) neuron model which was originally proposed to ac- 
count for the property of squid giant axons [12j and has 
been generalized with modifications of ion conductances 
[T^|. It describes the spiking behavior and refractory 
properties of real neurons and serves as a paradigm for 
spiking neurons based on nonlinear conductance of ion 
channels 01 • The dynamics of HH model is described 
by the following set of coupling ordinary differential equa- 
tions: one for the membrane potential V and the other 
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where £(t) is Gaussian white noise, D and Td is the in- 
tensity and the correlation time of the colored noise, re- 
spectively. In the following numerical simulations we 
took Td = 2msec. Differential equations given by Eqs. 
(|lafl - (jTe|) were integrated by the forth-order Runge-Kutta 
method. The colored noise in Eq. was solved follow- 
ing ideas by Fox et al. 0|. The time step was At = 0.02 
msec. 

We use spike-train rather than membrane potential to 
represent the temporal sequence of firing events of neu- 
rons. Because it is believed that most of information in 
neural systems is coded in the time sequence of action po- 
tential 17]. The spike-train is a binary time series with 
a value 1 at the time of action potential generations and 
at other times. The spikes were defined as 20mV-level 
crossing from below to above values of the membrane 
potential V. 
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FIG. 1: The waveform of the input periodic signal and three 
simulated spike-trains from the HH neuron, (a) The waveform 
of the signal used in simulations. The coordinates of signal 
was ignored, (b) The spike-train obtained by inputting strong 
signal (A\ — 2.0fiA/crn 2 ) into the HH neuron and noise ab- 
sent, (c) The spike-train produced by added noise (D = 5) 
to the strong sine signal used to obtain (b). (d) The spike- 
train recorded from simulations using the sum of the weak 
sinusoidal signal (Ai = 1.0fiA/cm 2 ) and noise (D = 5.0) as 
stimulus. 
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FIG. 2: (a) The waveform of the signal which is the sum 
of two sine wave which have frequency fi = 0.05kHz and 
/2 = 0.065kHz, respectively, (b) The spike-train produced by 
the HH neuron which received the combined signal without 
noise, (c) The spike-train recorded from the stochastic HH 
neuron (D = 10.0) responding to the combined signal (A2 = 
l.0fiA/cm 2 ). 

Some examples of spike-trains produced by the com- 
puter simulations of HH neurons responding to periodic 
signal with or without noise were presented. First, we 
considered a pure tone sinusoidal external stimulus: 

Jexi = Aisin(2nft). (3) 

In Fig. n we plotted the waveform of the pure tone sig- 
nal with frequency / = 0.05kHz followed by three spike- 
trains recorded from three trials with different conditions. 



The curve in Fig. Qfa) represented precisely the phase 
of sinusoidal signals used in different trials, but inten- 
sity was not included. When the sine signal is strong, 
Ai = 2.0/iA/cm 2 , and noise is absent, the spike-trains in 
Fig. njb) was activated. There is a spike in each period 
of signal. And the spikes fire precisely at a specific phase 
of sine signal during every positive half-cycle. There is 
a common delay between spikes and corresponding wave 
crests, which does not affect the temporal structure of the 
spike-train and the correlation between the input signal 
and the spike-train. According to the phase locking of 
firing to sinusoidal signal, we simplified the description 
of the correlation between the input signal and the spike- 
train as that every wave crest of the pure tone sine wave 
activates a spike. In the second trial, we added noise into 
the HH neuron which received the sinusoidal signal used 
upper. The intensity of noisy current was D — 5.0. The 
spike-train recorded from the simulation was shown in 
Fig. njc). Notice that most of the spikes deviated from 
the positions contrasting to that in the above spike-train 
and few spikes disappeared randomly. The size of devia- 
tion is random and small. In the third trial, we decrease 
the amplitude of periodic signal to A% — 1.0/iA/cm 2 , 
and remained the intensity of noise as D = 5.0. The 
spike-train computed from the trial was illustrated by 
Fig. IHd). Under this condition, the number of spikes 
had a sharp decrease compared with the spike-train ac- 
tivated by strong signal shown in Fig. [Jc). It is ran- 
dom that whether a wave crest activates a spike. The 
spikes in the spike-train also deviate from the positions 
of the wave crests. Although randomness entered the 
spike-train when noise was presence, the phase locking 
was affected slightly by noise. 

Pure tone stimulations have been used in many studies 
of neuronal response |g . Subsequently, we considered the 
signal which is a sum of two sinusoidal functions with 
different frequencies: 

Iext = A 2 {sin(2nfit) + sin{2nf 2 t)). (4) 

We took the frequencies of two sinusoidal functions as 
/1 = 0.05kHz, f 2 = 0.065kHz. In Fig. 0the waveform 
of the combined signal and the spike-trains activated by 
it were illustrated. In Fig. 0{b), during a period of the 
input signal, from 600 msec to 800 msec, three spikes 
fired and locked to the three highest wave crests. When 
noise was added to the neuron, the number of spikes in- 
creased. Some spikes fired randomly at the small neigh- 
borhoods of weaker wave crests. The stronger a posi- 
tive half-wave is, the larger probability it gives rise to a 
spike. According to the correlations between the input 
signal and the spike-train, we propose an interpretation 
of the encoding mechanism. The phase locking phenom- 
ena suggest that the response to periodic signal can be 
considered as a waveform sample: when a sensory neuron 
receives the periodic signal, it samples the input signal 
through the nonlinear response and the signal was sam- 
pled as the action potential. The sampling method has 
the following properties, (i) Sample points appear in the 
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FIG. 3: (a) The power spectrum of the spike-train excited 
by the weak pure sinusoidal current (A = 1.0fiA/cm 2 ) with 
noise (D = 5.0). (b) The power spectrum of the reconstructed 
spike-train produced by the MC simulation of the waveform 
sample method with I c = 2.0fj,A/cm and a = 2.0msec. 



neighborhoods of the wave crests of input signals, (ii) 
It is unnecessary that every wave crest is sampled. The 
probability that a wave crest is sampled depends on the 
strength of it. (iii) The positions of the sampling points 
distribute randomly in the neighborhoods of the wave 
crest of signals, (iv) The values of the sampling points 
do not depend on the signals which are sampled. All 
samples have the same value. 

To demonstrate the waveform sample method, we used 
it to recover the basic characteristics of neural signals. 
The power spectra of the spike-train was used to evalu- 
ate the frequency information of neural signals and assess 
the coherence of the spiking activity with the signal fre- 
quency, which has been used extensively in the studies of 
stochastic resonance as mentioned above. We computed 
the power spectral density from simulating the process 
that an individual HH neuron responded to periodic sig- 
nals and the Monte-Carlo (MC) simulations of the wave- 
form sample method. The power spectra were computed 
from 200 averages of the power spectra for the time series 
of the spike-train with length 2 17 using the fast Fourier 
transform [T^ |. 

Let us first consider the scheme that a stochastic HH 
neuron receives the pure tone sinusoidal signal described 
by Eq. (0 ■ We took the amplitude and frequency of sig- 
nal as Ay = 1.0^A/cm 2 and / = 0.050 kHz. This signal 
is too weak to activate HH neuron without noise (the crit- 
ical amplitude is A c — 1.52/j,A/cm 2 when / = 0.050/cHz). 
We calculated the power spectrum of the spike-train in 
the presence of noise {D = 5.0). The power spectrum 
was shown in Fig. |2{a). In the power spectrum, the main 
feature due to the periodic signal is the narrow peak at 
the fundamental frequency riding on a broad noise back- 
ground. Additionally, the harmonic of the fundamen- 



tal frequency exists in the power spectrum. This result 
agrees excellently with the experiment in Ref. PJ. The 
power spectrum quantify the information content of the 
spike-train without invoking any specific encoding mech- 
anism. The peak of power spectral density located at the 
fundamental frequency proves that the information of the 
period of the signal must be encoded in the spike-train 
by the individual HH sensory neuron. As an example, 
we consider a simple rule of the waveform sample. The 
probability that a wave crest gives a sample point is pro- 
portional to the amplitude of signal: 



Pi = h/Ic, 



(5) 



where, Ii = Ax is the amplitude of the pure tone sinu- 
soidal signal, I c is constant which is larger than Ax- The 
positions of sample points were subjected to Gaussion 
distribution: 



P(Si) 



1 



27TCT 



exp(- 



2a 2 



), 



(6) 



where U is the position of the zth wave crest, Si is the 
position of sample point that locks in the neighborhood of 
ti. We made MC simulations using the rule described by 
Eqs. © and © with I c = 2.0^A/cm 2 and a = 2.0msec. 
The power spectrum obtained from the MC simulation 
was compared with that of the HH neuron in Fig. [3] 
Although the background noise in the power spectra are 
different, the frequency structure of HH neuron has been 
recovered virtually. The difference between background 
noise in Fig. &b) and that in Fig. |2Ia) is due to the 
simple rule of waveform sample which does not involve 
the dynamics of HH neuron. 

Then we consider the combined signal described by Eq. 
QJ. The parameters were taken as A2 = 1.0/j.A/cm 2 , 
fx = 0.05kHz, and f% — 0.065kHz in simulations. Un- 
der this condition, the HH neuron can be activated by 
the combined current without the assistance of noise. In 
the power spectrum, as illustrated by Fig. 01a), the fun- 
damental frequencies give two high and narrow peaks 
which were marked by the labels using their positions 
in frequency domain. Also lots of additional frequency 
components arise and some of which are as strong as fun- 
damental frequencies. The frequency structure exhibited 
in this power spectrum was given by the precisely phase 
locking spike-train shown above in Fig. OJb). When noise 
was present, we simulated the process that HH neuron re- 
sponds to the combined signal. In the case of noise inten- 
sity was taken as D = 10.0, the power spectrum was com- 
puted and shown in Fig. Elb). The frequency structure 
has a observably change compared with Fig. 0fa). The 
fundamental frequencies fx and fi give the highest peaks 
in the power spectrum. Most of the additional frequency 
components disappear. The surviving additional compo- 
nents locate at fi — fx, 2/i and /1+./2 with weak strength. 
Noise affects the frequency structure of the output spike- 
train of sensory HH neuron and plays an optimized role in 
distinguishing the fundamental frequency components in 
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FIG. 4: (a) The power spectrum recorded from the simu- 
lation of HH neuron responding to the periodic signal de- 
scribed by Eq. Q without noise. The parameter were taken 
as fx = 0.05kHz, f 2 = 0.065kHz and A 2 = 1.0fiA/cm 2 . 
(b) The power spectrum produced by stochastic HH neuron 
(D = 10.0) responding to the combined signal, (c) The power 
spectrum obtained from the MC simulation of waveform sam- 
ple. 

the power spectrum. In the MC simulation of the wave- 
form sample of the combined periodic signal, the sample 
rule described by Eqs. J5j and © was used again. The 
symbol Ii was the height of the ith wave crest of the 
combined periodic signal. We took the parameters as 



I c = 2.0fiA/cm 2 and a — 3.0msec. The power spectrum 
of the waveform sample was shown in Fig. Efc), and it 
got the primary characters of the frequency structure of 
the Fig. E{b). The result demonstrates that the ran- 
domness in the waveform sample method captures the 
effect of noise in sensory neuron. The sample rule used 
in the MC simulation was too simple and crude to depict 
the neural process exactly, but the main characteristics 
of spike-trains of the HH model have been obtained from 
the waveform sample. The waveform sample method cap- 
tures the nature of the encoding mechanism that sensory 
neuron responds to periodic signals. 

The present work is the further investigation of the 
stochastic process theory of firing events (3, where the 
frequency structure of power spectra were described ele- 
gantly. We obtained the subtler correlation between the 
spike-train and the input signal, and proposed a wave- 
form sample method as an intuitional interpretation of 
the encoding mechanism. The waveform sample method 
is different from the conventional sample method used 
for processing signals, such as the even interval sample 
used for discrete Fourier transform |l9j | which does not 
give rise to the additional frequency components in power 
spectra. If the frequency of the signal received by neuron 
is low, the neuron will fire more than once in a period of 
signal. On the other hand, the dynamics of HH neuron 
can not lock the spike to every wave crest of the signal in 
the high frequency region because of the refractory prop- 
erties. Although it is not a proper encoding method of 
neurons in low and high frequency regions, it can also 
convey the frequency information into the sample points 
in all these frequency regions. Additionally, the waveform 
sample method suggests a method for reconstructing the 
spike-train and proposes a chance to utilize the neural 
encoded mode without invoking neuronal models. 
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